function error_probability()
        
    N = double([128, 64]);
    t = double(2);
    
    disp(rs_code_p_out(double(1e-7), double(64), double(2)));
    
    p_in = double((1:0.01:100)*1e-9);
    p_out = double(zeros(length(N), length(p_in)));
    
    for j=1:length(N)
        for i = 1:length(p_in)
            % p_out(j, i) = rs_code_p_out(p_in(i), N(j), t);
        end
    end
           
   %  plot(p_in, p_out(1,:), p_in, p_out(2,:));
    
end

function Perr = rs_code_p_out(pin, N, t)
    % Вероятность ошибочного 10битного кодового слова -> 
    %       врероятность получения ошибочного бита
    Pb = 1 - (1 - pin)^10;
    
    sum = 0;
    for i=0:t
        k = (nchoosek(N,i)*(Pb^i))*((1-Pb)^(N-i));
        sum = sum + k;
    end

    Perr = 1 - sum;
end